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Abstract 



We study the direct detection rate for SUSY cold dark matter (CDM) pre- 
dicted by the minimal super symmetric standard model with universal bound- 
ary conditions and large values for tan(3. The relic abundance of the lightest 
supersymmetric particle (LSP), assumed to be approximately a bino, is ob- 
tained by including its coannihilations with the next-to- lightest supersymmet- 
ric particle (NLSP), which is the lightest s-tau. The cosmological constraint 
on this quantity severely limits the allowed SUSY parameter space, especially 
in the case the CP-even Higgs has mass of around 114 GeV . We find that for 
large tan/3 it is possible to find a subsection of the allowed parameter space, 
which yields detectable rates in the currently planned experiments. 
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In recent years the consideration of exotic dark matter has become necessary in order 
to close the Universe JT]. Furthermore in in order to understand the large scale structure 
of the universe it has become necessary to consider matter made up of particles which were 
non-relativistic at the time of freeze out. This is the cold dark matter component (CDM). 
The COBE data || suggest that CDM is at least 60%. On the other hand during the 
last few years evidence has appeared from two different teams, the High-z Supernova Search 
Team || and the Supernova Cosmology Project Q which suggests that the Universe may 
be dominated by the cosmological constant A. As a matter of fact recent data the situation 
can be adequately described by a baryonic component Qb = 0.1 along with the exotic 
components flcDM = 0.3 and fl^ = 0.6. 

In fact the DAMA experiment has claimed the observation of one signal in direct detection 
of a WIMP, which with better statistics has subsequently been interpreted as a modulation 
signal || . 

The above developments are in line with particle physics considerations. Thus, in the 
currently favored supersymmetric (SUSY) extensions of the standard model, the most natural 
WIMP and candidate for CDM is the LSP, i.e. the lightest supersymmetric particle. In the 
most favored scenarios the LSP is the lightest neutralino, which can be simply described as 
a Majorana fermion, a linear combination of the neutral components of the gauginos and 
Higgsinos flj], |§- [jOD . Its stability is guaranteed by imposing R-parity conservation , which 
implies that the LSP's can disappear only by annihilating in pairs. Since this particle is 
expected to be very massive, > 30GeV, and extremely non relativistic with average 
kinetic energy T < lOOKeV, it can be directly detected flOfl - |TT| mainly via the recoiling 
of a nucleus (A,Z) in the elastic scattering process: x + (A ~^ X + (A %)* 

In order to compute the event rate one needs the following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level obtained in the frame- 
work of supersymmetry as described , e.g., in Refs. and Bottino et al ||. 

2) A procedure in going from the quark to the nucleon level, i.e. a quark model for the 
nucleon. The results depend crucially on the content of the nucleon in quarks other than u 
and d. This is particularly true for the scalar couplings as well as the isoscalar axial coupling 

0-0- 

3) Compute the relevant nuclear matrix elements, using as reliable as possible many body 
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nuclear wave functions. A complete list of references can be found in recent publications 



In the present work we are going to be primarily concerned with the first item of the 
above list. We will be concerned, in particular, with SUSY models predicting a LSP with a 
relic abundance of cosmo logical interest. Following the considerations of on the com- 
position of the energy density of the universe in scenarios with vanishing and non vanishing 
cosmological constant we assume flisph 2 in the range: 

0.09 < n LS ph 2 < 0.22 (1) 

Combining items 1) and 2) we can obtain the the LSP-nucleon cross section. We will 
limit ourselves in the experimentally interesting range: 

4 x 10~ 7 pb < ff™t m < 2 x 10~ 5 pb (2) 

We should remember, however, that the event rate does not depend only on the nucleon 
cross section, but on other parameters also, mainly on the LSP mass and the nucleus used 
in target. 

The calculation of this cross section has become pretty standard. One starts with rep- 
resentative input in the restricted SUSY parameter space as described in the literature, e.g. 
Bottino et al. ||, and Arnowitt et al. |9|], |TL], and recently f22]| . 

In the present work will show how the more restrictive version of the MSSM, with minimal 
supergravity (mSUGRA) and gauge unification, can predict an LSP consistent with a relic 
abundance in the range ([]]) and experimentally interesting detection rates. The space of 
SUSY parameters compatible with a CP— even Higgs mass nth ~ 114 GeV will restrict the 
cosmological preferred area to the region where \ — f coannihilations effects are sizeable. 



We will consider the MSSM which described in detail in Ref. |2(J and closely follow the 
notation as well as the renormalisation group (RG) and radiative electroweak symmetry 
breaking analysis of this reference. The full one-loop radiative corrections to the effective 
potential for the electroweak symmetry breaking and one-loop corrections to certain particle 
masses which have been evaluated in Ref. [^3[] (Appendix E) are incorporated as described 



in Ref. |24| . We do not impose any unification condition to the Yukawa couplings, however, 



even though we maintain our calculation in the large tan(3 regime. We incorporate in our 



2 



calculation the two-loop corrections to the CP-even neutral Higgs boson mass matrix by 
using the FeynHiggs subroutines Ref. [^3J this improvement is significant, especially for the 
lighter eigenvalue, at large values of tan/3. 

We restrict our a study to the case of universal soft SUSY breaking terms at the Grand 
Unification scale, Mqut, i-e., a common mass for all scalar fields mo, a common gaugino 
mass M1/2 and a common trilinear scalar coupling A . Our effective theory below Mqut 
then depends on the parameters (/i = fi(M GUT )) 

m , M1/2, A*o, A , «g, M GUT , h t , , h b , h T , tan/3 , 

where a G = 9gI^ (9g being the GUT gauge coupling constant) and h t ,hb,h T are respec- 
tively the top, bottom and tau Yukawa coupling constants at Mqvt- The values of a G and 



Mqut and |/xo| are obtained as described in Ref. |20|j . 

For a specified value of tan (3 at the common SUSY threshold, defined as M s = ^/m^m^, 
we determine ht at Mqut by fixing the top quark mass at the center of its experimental range, 
Wt(m^) = 166 GeV. The value of h T at Mqvt is fixed by using the running tau lepton mass at 
mz, m T (mz) = 1.746 GeV. We also incorporate the SUSY threshold correction to m T (Ms) 
from the approximate formula of Ref. |23[ . 

The value of hb at Mqut is obtained such that 

m b(mz)Wf = 2.90 ± .14 GeV (3) 

in agreement with the experimental value given in JZE]. The SUSY correction |27| to the 
bottom quark mass is known to increase with tan f3. This correction originates mainly from 
squark/gluino and squark/chargino loops and has the same sign as /i (in our convention, 
which is the oposite to the one used in Ref. p0| , ^3l^4l ). With the values of the SUSY 
parameters presented in Fig. 1, the SUSY correction to rrib decreases by a 4-5% as goes 
from the lower to the higher value in the considered range. Following the analysis of ||28|| , 
the values of eq. ^| correspond to mb{™b)sM — 4.2 ± .2 GeV which is in the range quoted in 
Ref. |9|. 

The choice \x > leads to a constraint on the parameter space coming from the lower 
bound on b — * sj. As a result the relatively light vlaues for and the obtained detection 
rates are suppressed. From the experimental results j3(J we consider a range within the 95% 
C.L. : 
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2.0 x 1(T 4 < BR(b -> 57) < 4.5 x 10" 4 (4) 
Our determination of BR(b — > 37) follows the procedure described in ref. P4" |. However 



appropiate calculations for the next-to-leading order (NLO) QCD corrections to the SUSY 
contribution at large values of tan/3 have become available |]3T|. The SUSY contribution 
decreases as Ms increases, we found that the lower bound in (|4j) determines a lower limit 
for my.. Therefore in order to perform a consistent prediction for BR(b — > 57) we have 
implemented the NLO QCD correction from the first Ref. [|3T| to the analysis described 



in 



The GUT values of M1/2 and mo are traded by the value of the mass of the LSP and the 
mass splitting between the LSP and the NLSP (next to lightest supersymmetric particle), 
which in our approach is one of the f's. The parameter Af 2 = {m^ 2 — m^jm^ is one of our 
basic independent parameters. This is motivated by the role played by coannihilations x ~ t 
in the determination of the cosmologically preferred area of the parameter space. We keep 
Aq = in all our calculations. 

The composition of the LSP on the model under consideration can be written in the basis 
of the gauge and Higgs bosons superpartners as: 

X = X° = C n B + C 12 W + C 13 Ht + C U H 2 . (5) 

In the parameter space we study, x is mostly a gaugino with P = \Cn\ 2 + \Cu\ 2 > .9, with 
the Bino component being the most dominant one. 

The parameter space considered relevant for our analysis is displayed in Fig.|l]. For 
convenience we plot the Ms, vtia and mo versus the the LSP mass. Since in our analysis 
the LSP is mostly a bino, one finds My 2 ~ 2.5m^. Two values of tan/3 were considered, 
one higher and one lower, which we view as representative. The higher one, tan/3 = 52, 
corresponds approximately to the unification of the tau and top Yukawa couplings at Mqut- 
We should clarify, however, that in making this choice we do not adhere to any particular 
Yukawa unification model. We find that the parameter space with this value of tan/? is 
representative for the purpose of this work . The lower value, tanf3 = 40, results in a cr" c ^f a e r on 
smaller by one order of magnitude as we will see later. The shaded areas correspond to the 
range of values taken by mi, Ms as ranges from to 1. The area associated with mo 
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for the same range of A^ 2 is wider as shown by the dashed and solid lines. Lower values of 
tan/3 will tend to increase the values for tua- This will lead to a decrease of the detection 
rates as we will see later. 



It can be shown iTTH, [181, iT3fl that the non directional rate can be cast in the form: 



R = R t(a, Qmin)] [1 + h(a, Q min )cosa)] (6) 

where a = \jjL r (A) b _1 with b the nuclear size parameter, fj, r (A) is the LSP-nucleus 

reduced mass and the parameter characterizing the LSP velocity distribution assumed here 
to be Maxwell-Boltzmann. The quantity h describes the modulation, i.e. the dependence of 
the rate on the Earth's motion, which is of no interest to us in the present work. Q m in is the 
energy transfer cutoff imposed by the detector. The parameter t takes care of the folding 
with the LSP velocity distribution and the structure of the nucleus and it also depends on the 
LSP mass. The most important to us parameter is R, which carries the essential dependence 
on the SUSY parameters, and it can be cast in the form: 

- p(0) m rr^rf, . v (u 2 )^ , 

R = ^A^y {v)[ ^ + + (7) 

Where p(0) is the LSP density in our vicinity, which will be taken to be 0.3 GeV/cm 3 , 
and m is the detector mass. Ej,z = S,V,spin are associated with the scalar, vector and 
spin contributions. In the present paper we will consider only the scalar contribution, which 
offers the best chance for measurable detection rates, in particular for medium-heavy nuclear 
targets, like 127 1. In this case 

S S = a £{A) {A 2 l(f s - fl^^f] ^ < x oA\^) 2 (8) 



where is the nucleon mass, p r (iV) is the LSP-nucleon reduced mass, is the scalar 
LSP-proton cross-section and 

a = _L(G FmiV ) 2 ~ 0.77 x KT 38 cm 2 (9) 

The coherent scattering can be mediated via the the neutral intermediate Higgs particles 
(h and H), which survive as physical particles. It can also be mediated via s-quarks, via the 
mixing of the isodoublet and isosinlet s-quarks of the same charge. In our model we find 
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that the Higgs contribution becomes dominant and, as a matter of fact, the heavy Higgs H is 
more important (the Higgs particle A couples in a pseudoscalar way, which does not lead to 
coherence) . As it has been pointed out in Ref . |2TJ large values of tan f3 and low values for the 
mass of the CP-even Higgs masses (m^ and ran) enhance the Higgs exchange amplitude. This 
is of particular interest in our analysis, since, as tan /3 increases, the electroweak symmetry 
breaking imposes lower values for for the pseudoscalar Higgs mass tha (see Fig. 1). This 
implies a lower value for of m#. The changes on are not so important, since its value can 
only move below an upper bound of about 120-130 GeV (see for e.g |25|). Since the values 



predicted for the /i-parameter also decrease, the Higgsino component of the LSP enhances 
the values of cr^f a e ° n . For all the values of the m-j, however, reported in the present work the 
condition P > .9 is maintained. It becomes even more stringent, P > .95, for > 100 GeV. 



It is well known that all quark flavors contribute, see e.g. Drees et al ||12|| , since the 
relevant couplings are proportional to the quark masses. One encounters in the nucleon not 
only the usual sea quarks (««, dd and ss) but the heavier quarks c, b, t which couple to the 
nucleon via two gluon exchange, see e.g. Drees et al and references therein. 

As a result one obtains an effective scalar Higgs-nucleon coupling by using effective quark 
masses as follows 

m u -> f u m N , m d -> f d m N . m s -> f s m N 
m Q ~~ > fQ m N, (heavy quarks c,b,t) 

where is the nucleon mass. The isovector contribution is now negligible. The parameters 
f q , q = u,d,s can be obtained by chiral symmetry breaking terms in relation to phase shift 



and dispersion analysis. Following Cheng and Cheng [|14]] we obtain 

f u = 0.021, f d = 0.037, f 8 = 0.140 (model B) 

/„ = 0.023, f d = 0.034, f s = 0.400 (model C) 

We see that in both models the s-quark is dominant. Then to leading order via quark loops 
and gluon exchange with the nucleon one finds: 

f Q = 2/27(1 - E q f q ) 
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This yields: 



f Q = 0.060 (model B), f Q = 0.040 (model C) 

There is a correction to the above parameters coming from loops involving s-quarks [|HJ and 
due to QCD effects. Thus for large tan/3 we find : 

f c = 0.060 x f .068 = 0.064, f t = 0.060 x 2.048 = 0.123, f b = 0.060 x 1.174 = 0.070 

(model B) 

f c = 0.040 x 1.068 = 0.043, f t = 0.040 x 2.048 = 0.082, f b = 0.040 x 1.174 = 0.047 

(model B) 



For a more detailed discussion we refer the reader to Refs [p~2| , p~3|1 . With the above ingredients 



and employing the obtained constraints on the SUSY parameter space one can proceed to 
calculate first the parameter fg associated with the scalar contribution. Then combine it 
with the effective Higgs-quark couplings f q and /q to obtain the nucleon cross section (the 
contribution to the scalar coupling arising from the L-R s-quark mixing is in our model 
negligible). 

The values for <T™ c ^f a e r on for the two nucleon models described above and the input SUSY 
parameters presented in Fig.l are shown in Fig. [2]. As we can see nucleon model C leads 
to higher cross sections. For tan /3 = 52 the effect of changing A ?2 (equivalently mo in the 
range shown in Fig.l) from to 1 has a small effect on the cross section (not sizeable in 
the figure) while for tan (3 = 40 the lines are widened and become bands. This effect can be 
understood in the light of Fig. 1. As tan/5 decreases the values of (and hence m# and 
the /^-parameter) increase and span a larger area as A ?2 varies from to 1. This leads to a 
substantial decrease in the cross section for the lower value of tan (3. 

The fact that (x) is mostly a B implies that the main contribution to its annihilation 
cross section arises from s-fermion (squark, s-lepton) exchange in the t- and u-channel leading 
to // final states (/ is a quark or lepton). If, however, the mass of % is close to the one 
of the NLSP, coannihilations between the two particles must be taken into account [J32| . 
The inclusion these coannihilation effects results in a dramatic reduction of the (x) relic 
abundance as the two lightest SUSY particles approach in mass [50,53151] . We estimate the 
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relic abundance of the LSP (x), by employing the analysis of Ref. 0] which is appropriate 
for large tan f3 and includes coannihilations \ — f , suitable for Bino like LSP. To the list of 



coannihilation channels of presented in |2C ] , one should add the t-channed sneutrino exchange 
in the W + W~ and H + H~ production. 

We should, at this point, clarify that in the parameter space considered here no resonances 
in the s-channels were found. In other words the s-channel exchange of A , h, H, Z into 
T2T2* never becomes resonant in the parameter space of our analysis. We can see in Fig. 1, 
however, that a line mass = 2m^ will be above of the region for the case of tan j3 = 52, 
while for tan/5 = 40 it will not. However we should emphasize here, that the position the 
rriA band displayed is Fig.l respect a line of mass = 2m^ is very sensitive to small changes in 
tan j3 and the values mt,rrib and the GUT values for Aq and m^. Therefore, at the large values 
of tan (3 it is possible to find sectors of the space of parameters where m^ m m^., in these cases 
the the adequate treatment of the Higgs mediated anihilation channels will be determining 
for an accurate calculation of Qlsp h 2 . The choice of parameter space in the two examples 
we present is aimed to illustrate the decisive role of tan j3 in the LSP detection rates. In the 
two scenarios we choose, anihilation resonant channels are not present and coannihilations 
are required in order to predict a cosmologically desirable LSP relic abundance. 

Fig. |3| shows the cosmologically allowed area from CDM considerations corresponding to 
the parameter space described in Fig.l. The b — ► 57 constraint will exclude the parameter 
space marked on the left of the figures. The larger range on the the values of the LSP 
masses corresponds to the narrow band of the mass splitting between the LSP-NLSP, which 
enhances the effect of coannihilations. We find the lower of the bounds of Eq. (JJ) in the 
displayed area corresponds Af 2 = .045 for tan/5 = 52 and Af 2 = .03 for tan/? = 40, the 
effect of the x ~ t coannihilation disappears for Af 2 rs .25. Higher values for m^ 2 , enhance 
the neutralino relic abundance. We found that for values of m^ > 105GeV, a choice of 
parameters such that m^ 2 > 2m^ leads to values for Qisph 2 > .22, the upper limit we 
consider as cosmologically desirable. 

The b — > 57 constraint is more restictive than the bounds on the masses of the f and 
chargino. Lines of constant mass for lighter f can be easily deduced from the coordinates of 
the graphs in Fig. |3], the ones corresponding to values of about 95 GeV will start intersecting 
the cosmologically favored area in both examples, they would lay in the area excluded by 
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the b — > s'y constraint. Lines of constant chargino mass would appear in the graphs of 
Fig. ^ as almost vertical lines. The line corresponding to m^+ = 110 GeV will be located 
at rriy, = 64 GeV for tan/3 = 52 and = 62 GeV for tan/3 = 40. In the space of 
parameters under consideration, the difference m^+ — is large enough to consider any 
effect of x + ~ X coannihilations in the x renc abundance. Other scenarios allowing sizeable 
X + ~ X coannihilation effects may be possible for large tan/3. Such a study, however, is 
beyond the scope of the present work. 

We find that the possible values of lightest neutral CP even Higgs mass, imposes a severe 
constraint on the cosmo logically allowed area. We also found that, in a fashion similar to 
the results presented in Ref. |35j for the low tan/3 scenario, in the large tan/3 case a value 
of rrih = 114 GeV will push the cosmologically allowed area to the values of the parameters 
where x ~ t coannihilations are important (see Fig. |3|). Changes in the values of Af 2 will 
not affect significatively the values of rrih since, as we can see in the shaded areas of Fig. 
[1], the changes in Af 2 imply small changes on the values of rriA and Ms- Conversely small 
changes on rrih are translated in significative displacements of the line of constant rrih shown 
in Fig. |3]. Indeed this line would be better understood as a band of approximately 20 GeV 
width on for changes on rrih from 113 GeV to 115 GeV (these changes can be understood 
from the shape of the curves shown, e.g, on the first reference of Ref. P5|]). We can note, by 
observing Fig. |l|, that under the conditions imposed on the parameter space larger values 
for tan/3 imply lower values for itia at the same range of m^. This helps to understand the 
different position of the lines of constant rrih in the two graphs of Fig. |3|. 

As we can see in Fig. |3], a value of rrih = 114 GeV will exclude the the space of parameters 
suitable for detection for tan/3 = 40 while for tan/3 = 52, the line corresponding to the 
same value for rrih lies inside the range of Eq. (0) for model C. The lines corresponding to 
rrih = 105 GeV would be also almost vertical lines located at the abscissa rriy. « 63GeV for 
tan/? = 52 and below the displayed area for tan/3 = 40. Therefore the dependence of our 
conclusions on rrih is clear. 

With the above ingredients one can proceed, as outlined above, with the calculation of 
the event rates. For a typical nucleon cross section, taken, e.g., from the set associated with 
tan (3 = 52, exhibited in Fig. ^|, we plot in Fig. ^ the event rate as a function of the LSP 
mass in the case of the popular target 127 1. Since the detector cut off energy is not so well 
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known, we consider two possibilities, namely Q m in — and 10 KeV. The values of the 
quantity t, see Eq. (|5|), employed here have been taken from our earlier work [19]. In these 
plots we have excluded LSP masses less than 130 GeV to avoid a logaritmic scale. They can 
easily be extrapolated to values of vn^sp down to 100 GeV. Anyway event rates in this LSP 
mass region are very large in our model and it unlikely that they would have been missed 
by the existing experiments. 

In summary, we have found that the most popular version of the MSSM with gauge 
unification and universal boundary conditions at the GUT scale, and a parameter space 
determined by large values of tan/5, can accommodate a cosmologically suitable LSP relic 
abundance and predict detection rates, which can be tested in current or projected experi- 
ments. 

We should mention, however, that the calculated detection rates can vary by orders of 
magnitude, depending on the yet unknown LSP mass. We will not be concerned here with 
uncertainties, which are essentially of experimental nature, like the threshold energy cut off 
imposed by the detector. This reduces the rates, especially for small LSP masses. In the 
LSP mass range considered here the reduction is about a factor of two (see Fig. |]). The 
next uncertainty comes from estimating the heavy quark contribution in the nucleon cross 
section. This seems to be under control. We take the difference between the models B and 
C discussed above as an indication of such uncertainties. They seem to imply uncertainties 
no more than factors of two (see Fig. §). The nuclear uncertainties for the coherent process 
are even smaller. 

We believe, therefore, that, concerning the direct LSP detection event rates the main 
uncertainties come from the fact that the SUSY parameter space is not yet sufficiently 
constrained. The parameter space may be sharpened by the accelerator experiments, even 
if the LSP is not found. We should mention here, in particular, the Higgs searches, since, 
as we have seen, the role of the Higgs particles in direct SUSY dark matter detection is 
crucial. It is not an exaggeration to say that the underground and accelarator experiments 
are complementary and should achieve a symbiosis. 

Note added: After the first version of this work was completed the Muon (g — 2) Col- 
laboration had published a new meassurement of the anomalous magnetic moment of the 
muon |3H||. The result presented differs from the SM prediction by 2.6a. If this discrepancy 
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is attributed to SUSY contributions, the (g — 2) may impose a further constraint on the 
parameter space presented here. When we consider the reported discrepancy at the 2 — a 
level as the allowed range for the SUSY contribution, the values we found do not have a 
significative impact in the paramter space considered in this work. Using the calcualtion 
presented in ref. |J7| the upper bound will exclude mostly the area already excluded by 
b — > S7, while the area excluded by the lower bound corresponds to values of > 300GeV. 
We are thus in agreement with Ref. |38j in the space of parameters we overlap. 

We would like to thank S. Khalil for useful discussions and for reading the manuscript 
and T. Falk for pointing out the missing coannihilation channels on the list of Ref. ||20|| . 
M.E.G. thanks the University of Ioannina for kind hospitality. This work was supported by 
the European Union under TMR contract No. ERBFMRX-CT96-0090 and TIENEA 95 of 
the Greek Secretariat for Research. 
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FIG. 1. The values of mo, and Ms as functions of m^ or for tan/3 = 40 (left) and 
tan/3 = 52 (right), /x > 0, Ao = 0, the upper boundaries (lower) on the shaded areas and the upper 
(lower) dashed line corresponds to m? 2 = 2 x (mf 2 = m^). 




m LSP (GeV) 

FIG. 2. Values for cr^^° n \ corresponding to the input parameters presented in Fig. [[]. The 
upper lines corresponds to tan/3 = 52, the solid (dashed) one was obtained using model C (B). The 
effect of changing with Af 2 from to 1 is not sizeable for tan j3 = 52. The bands are obtained using 
tan j3 = 40, the darker shaded band corresponds to model C with the upper bound corresponding 
to Af 2 = 1 and the lower to Af 2 = 0. The light shaded band corresponds to model B, for the same 
range of values of A^ 2 




FIG. 3. The cosmologically allowed region in the — Af 2 plane (shaded area) for tan (5 = 40 
(left) and tan/3 = 52 (right). The vertical bands in graph on the left correspond to the bound 
a ^smiwr° n ^ = ^ ' 10~ 7 p&, obtained from figure 2 for both models. The vertical lines on the graph of 
the right correspond to the bounds cr^™j^ on ^ = 4 • 10~ 7 p& (lines towards the right of the graph) 
and 2 • 10 _5 p6 for the models indicated. The marked areas on the left are excluded by b — > sj. 
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FIG. 4. The Total detection rate per (kg — target)yr vs the LSP mass in GeV in the case of 
127 1, corresponding to model B (thick lines) and Model C (fine lines). We used the parameter space 
corresponding to tan/3 = 52 and m/, > 114 GeV. On the solid curves we used no detector threshold 
energy cut off. On the dashed ones, the value 10 KeV was employed. 



